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Understanding the stability of thin film nanomagnets with perpendicular magnetic anisotropy 
(PMA) against thermally induced magnetization reversal is important when designing perpendic- 
ularly magnetized patterned media and magnetic random access memories. The leading-order de- 
pendence of magnetization reversal rates are governed by the energy barrier the system needs to 
surmount in order for reversal to proceed. In this paper we study the reversal dynamics of these 
systems and compute the relevant barriers using the string method of E, Vanden-Eijnden, and Ren. 
We find the reversal to be often spatially incoherent; that is, rather than the magnetization nipping 
as a rigid unit, reversal proceeds instead through a soliton-like domain wall sweeping through the 
system. We show that for square nanomagnetic elements the energy barrier increases with element 
size up to a critical length scale, beyond which the energy barrier is constant. For circular elements 
the energy barrier continues to increase indefinitely, albeit more slowly beyond a critical size. In 
both cases the energy barriers are smaller than those expected for coherent magnetization reversal. 



Thin film elements with perpendicular magnetic 
anisotropy (PMA) can have magnetization directions 
that are thermally stable at room temperature at the 
nanometer scale, a feature that makes them useful in 
information storage and processing, for example in pat- 
terned media [l[ and spin-transfer MRAM [2]. A key 
issue in determining stability is how the energy barri- 
ers and transition states of these elements depend on 
their lateral size. For elements larger than the exchange 
length (typically ~ 5 nm in transition metal ferromag- 
nets) the assumption of coherent reversal of the mag- 
netization breaks down and the transition state is not 
uniformly magnetized. Due to the multiscale character 
of micromagnetism, analytical calculations are compli- 
cated and transition states have been calculated only for 
a handful of physical systems Numerical calcula- 

tions are usually necessary for the majority of systems. 
In this paper, we use the string method [6j] to find the 
transition states and activation energies in thin film ele- 
ments with PMA. 

Transition state theory indicates that reversal occurs 
through states corresponding to critical points of the 
magnetization energy functional, where Vm^ = 0. The 
probability for thermally induced magnetization reversal 
is expected to follow the Arrhenius law e - u / k & T where 
U is the energy difference between the transition state 
and the met ast able configuration 0. Here, we iden- 
tify the minimum energy paths (MEPs) on the energy 
landscape, which allow us to determine the transition 
states; we then study the dependence of these states, 
along with their corresponding energy barriers, on sam- 
ple size and geometry. Our two main conclusions are: 
(i) typically, models that assume uniform magnetization 
seriously overestimate the magnitude of the activation 
energy, especially in larger systems; and (ii) even in sit- 
uations where the assumption of uniform magnetization 
is valid, that is, in smaller size systems, finite size effects 
have a strong impact on the value of the energy barrier. 



The micromagnetic energy of this system is given by 



[8]: 
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Here m is the normalized magnetization vector at po- 
sitions r and the integral is evaluated over the sample 
volume V = Clt, where Cl is the cross sectional area of 
the sample and t its thickness. The first term under 
the integral is the exchange energy, with exchange con- 
stant A; this term favors homogeneous magnetization in 
the sample. The second term is the magnetocrystalline 
anisotropy with constant K, which favors specific ori- 
entations of the magnetization. The third term corre- 
sponds to the Zeeman energy, which favors magnetiza- 
tion aligned with the external magnetic field H z : M s is 
the saturation magentization and /no is the permeability 
of the vacuum. The last term in Eq. [T] represents the 
demagnetizion energy; it captures the long-range dipole- 
dipole interactions between different regions of the sam- 
ple: 
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where m' is the magnetization at position r'. 

The macrospin approximation corresponds to A — >• oo, 
so that the sample magnetization rotates uniformly. In 
this case, m and m' can be factored out of the integrals. 
The dipole-dipole interaction then simplifies to E dd = 
^fioMgtQ m N m, where the demagnetizing tensor N 
is defined to be 
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For very thin and extended films (Cl —> oo, t = 0) the 
off-diagonal components of N tend to zero and the diag- 



onal terms approach N z 



1, N x 



0, N yy 0. In 
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this limit, the magnetostatic energy density simplifies to 
Edd — \\i^M s tVtm 2 z ^ and the in-plane components of the 
magnetization (m x ,m y ) are assumed to have negligible 
contribution to the magnetostatic energy. 

Summarizing, in the macrospin model for thin and ex- 
tended films, the magnetic energy reduces to E m , defined 
as: 
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where m z can have values between ±1 and the compo- 
nents m x and m y no longer play a role. An important 
conclusion of this work is that even in cases of very small 
aspect ratio tj y/Q this approximation is invalid, resulting 
in an underestimation the value of the energy barrier. 

To calculate the MEPs associated with the full mag- 
netization energy (pQ) and thereby determine the corre- 
sponding transition states and energy barriers, we uti- 
lize the String Method in conjunction with OOMMF (<|. 
We use 100 images between the two energy minima and 
reparametrize the string every 40 ps. Our initial guess 
path passed through a fully randomized magnetization 
configuration and was allowed to evolve according to the 
String Method prescription to its minimal energy path. 
Details are presented elsewhere (Tol-fl2j| . 

Using this approach, we studied circular and square 
layers of thickness t = 1.6 nm for a variety of diame- 
ters and side lengths L (the area Q is of order L 2 ). We 
consider a material with saturation magnetization M s = 
713 x 10 3 A/m, exchange constant A = 8.3 x 10" 12 J/m, 
and anisotropy constant K = 403 x 10 3 J/m 3 -parameters 
similar to those of Co/Ni thin film nanomagnets stud- 
ied experimentally in Refs. jl3l - [l5l |. A constant field 
perpendicular to the film fioH z was applied with a max- 
imum magnitude equal to the coercive field, defined by 
HqH k = (2K - fi M 2 )/M s = 0.245 T. Any finite ex- 
ternal field will break the degeneracy between the two 
lowest lying magnetization states (m z = ±1). Any field 
larger than Hk will deterministically switch the magneti- 
zation, owing to an induced instability of the higher-lying 
in energy metastable (for H < Hk) state. 

We investigate the transition from the metastable state 
(hereafter that with downward magnetization) to the 
ground state (upward magnetization) and compare with 
the results from the macrospin model in which we use 
N zz — 1 except when indicated otherwise. A typical re- 
sult is presented in Fig. [1] where the reversal can be seen 
to occur by propagation of a Bloch wall across the sam- 
ple. 

The field dependence of the activation energy is shown 
in Fig. [2l As expected, the activation energy decreases 
as the field approaches the coercive field. 

The activation energy of the macrospin model (U m ) 
can be obtained using Eq. (jlj. We first find the transition 
state by solving = ; we find that the transition 
occurs at m z = —H z /Hk- The difference in energies 
between the transition state m z = —H z /Hk and the 
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Figure 1: Images of the minimum energy path. Images shown 
are at locations 0, 20, 40, 60, 80, 100 in the string. Red and 
blue (light and dark) represent downward and upward mag- 
netization respectively. Reversal occurs by nucleation of a do- 
main on one corner, and propagation of a domain wall across 
the material. The transition state is located close to image 
20. The curve presents a kink near the 55th image. This 
occurs when the traveling domain wall reaches the opposite 
corners of the element and this feature has no effect on the 
magnitude of the transition rate. 
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Figure 2: Field dependence of activation energy for ele- 
ments with two different sizes. The smaller samples follow 
more closely the predictions of the usual macrospin model 
(N zz = 1). For large samples, it can be seen that the acti- 
vation energy is lower than that predicted by the macrospin 
model (N zz = 1) for a large range of values of the external 
field. In both cases the small value of A permit nonuniform 
magnetization configurations and the smallest sample behaves 
more as a macrospin. 
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This relation is shown in Fig|2] and compared to our nu- 
merical results. For large samples the activation energy is 
clearly lower than the prediction of the macrospin model, 
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Figure 3: Comparison of two versions of the macrospin model: 
N zz — 1 corresponds to the commonly used thin film limit and 
N zz — 0.95 is a film with finite aspect ratio. Blue circles are 
calculated with the string method for a very large exchange 
constant. 



as the macrospin does not account for spatial variation 
of the magnetization. 

Two additional points needs to be made regarding the 
way the energy barrier changes with the size of the ele- 
ment. First, we verified that the activation energy bar- 
rier is independent of the cell size used in the numerical 
integration scheme once the cells become smaller than 
the exchange length; this is a standard requirement of 
numerical micromagnetics. Therefore we conclude that 
our numerical results, such as those shown in Fig. 2, are 
not affected by numerical artifacts induced by changing 
the relative length scales of the system. Second, an ef- 
fective reduction on the element size can be achieved 
by increasing the exchange constant to a large value 
(A = 8.3 x 10 _10 J/m). In this limit our result should ap- 
proach the macrospin model even in our largest samples. 
The results are shown in Fig. [3j where the activation en- 
ergy is compared to two versions of the macrospin model 
with different magnetostatic tensor components. Our re- 
sults indicate that merely increasing A is not sufficient to 
bring the macrospin model into agreement with results 
found using the String Method. A further refinement of 
the magnetostatic energy term, in which the film is not 
treated as an infinite plane, is necessary. In this case, the 
demagnetization tensor is calculated using formulas due 
to Newell [16] for values t = 1.6 nm and L = 100 nm. 
Using these results yields a good agreement. This shows 
that in the macrospin model finite size effects are impor- 
tant even for very small aspect ratio (t/VU) thin film 
geometries. It must be emphasized that the gain in ac- 
curacy is obtained without a significant increase in com- 
putation time. 

The finite size effect in the limit of large A can be 
explained as follows: if the lateral sample dimensions 
are comparable with the domain wall width, there is a 
buildup of surface charge in the lateral faces of the sam- 
ple as the magnetization reverses; the sample is effec- 
tively a finite rectangular box. For intermediate values 
of A, such that the film can be considered to be wide 



in the plane but vertically thin (when compared to the 
exchange length), the domain wall width becomes much 
smaller than the lateral dimensions and the wall does 
not intersect with the sample sides and produce no sur- 
face charges. Consequently, there is no contribution to 
the magnetostatic energy and the sample is effectively 
an infinite plate. In the infinite plate limit, Gioia and 
James [17] nave argued that the magnetostatic energy 
can be absorbed into an effective crystalline anisotropy 
K' = (i^-^). 

Size effects on the activation energy in square ge- 
ometries are summarized in Fig. 2J where two regimes 
can clearly be identified. For small samples, the ac- 
tivation energy increases with the size of the element; 
here the macrospin model holds. Beyond a certain size 
(~ 50 nm), however, the activation energy becomes es- 
sentially constant. For large elements, the radius of 
the reversed nucleus is independent of the element size. 
For all sizes considered, the reverse domain resembles a 
circular region with domain walls making right angles 
with the sample's edge (see second image at Fig. [I]). 
At a given configuration Eq. (pQ) can be approximated 
as the sum of three dominant terms: a switched do- 
main with area = (ttR 2 /4) with energy density per 

unit area £i approximate to t ^°^ s — K — HQH z M^j] 
an unswitched domain with downward magnetization of 
area ^2 = ^ — ^1 = ^ — ttR 2 /4 and surface energy 
density £2 ~ t ^°^ s — K + /noH z M s ^j ; and the domain 
wall of length S3 = 2nR/A with a linear energy den- 
sity per unit length A3 « 4tvA^[8|. The total energy 
£1^1 + £2^2 + $3 A3 is maximum for R* = 39 nm, with 
an activation energy of U ~ 16 x 10 _20 J. This is a good 
approximation to the energy plateau presented in Fig. 2J 
Our estimate for R* predicts to good approximation the 
onset of the plateau region in the activation energy: the 
samples must be large enough to accommodate the nu- 
cleus of size R* at a distance of at least a few exchange 
lengths from the corners of the sample. A direct mea- 
surement from the numerical magnetization configura- 
tion gives R* = 37 nm. 

The w idth of the domain wall is expected to be 
2yJ 'A/ ' K' = 20 nm; we obtain a numerical value of 16 nm. 
We expect that the narrowing of the domain wall is the 
result of dipole-dipole interactions at the interface be- 
tween the two regions. Since the magnetostatic interac- 
tion favors antiparallel configurations, the spatial rate of 
twist perpendicular to the domain wall increases. The in- 
crease in exchange energy is compensated by a decrease 
of magnetostatic energy. 

The optimum radius R* can be obtained from the 
above prescription to be: 

R* =le^VQ^l (6) 



where Z ex = y j^w^ * s tne exchange length and Q — 
2 -j^2 is the uniaxial material parameter. The optimum 
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Figure 4: Activation energy for square samples (blue/dark) 
and circular samples (red/light). Inset: energy barriers for a 
square of side 100 nm for different values of K, the points are 
obtained with the string method and the line obeys Eq.©. 



radius depends on H z and diverges as H z —> 0. This is 
reasonable since in the absence of an external field there 
is no preferred magnetization direction, and so the en- 
ergy barrier is produced by a straight wall parallel to one 
of the sides of the sample. However, in order to obtain a 
reasonable estimate of the energy barrier we need to im- 
pose two conditions. First, we must guarantee that the 
switched domain fits in the sample by R* < L; second, 
we require the domain wall width to be smaller than the 
square side 2y/A/K' < L. Rewriting 2^/A/K' as V ex 



and using Eq. (|6]) we arrive at the following condition 
on the external magnetic field for our model to be valid: 

H z > 2M S (^ff ) • For example, with L = 50 nm, the con- 
dition is fioH z > 0.02 T. When this condition is satisfied 
the activation energy can be estimated by 
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The behavior of the energy barrier for square samples 
results from the presence of sharp corners that favor nu- 
cleation. In contrast, circular shapes do not present such 
behavior. Fig. 2] shows the computed activation energies 
for circular samples as a function of in plane diameter. 
The activation energy continues to increase as the diam- 
eter increases. As a result circular samples become more 
thermally stable for increasing size. 

We tested the validity of Eq. ([7j) by varying the value of 
the magnetocrystalline anisotropy. The results are repre- 
sented on the inset of FigHJ There is a clearly increasing 



and nearly linear behavior, as predicted by the analytical 
expression. 

In conclusion, we used the String Method to deter- 
mine activation energies for thermally induced magneti- 
zation reversal and compared the results to a macrospin 
model. Our results indicate that finite size effects on the 
macrospin model are important even for very thin sam- 
ples. The energy barrier of the macrospin model can be 
calculated more accurately and at no additional compu- 
tational cost by using the full demagnetization tensor. In 
addition, the activation energy for square samples as a 
function of size displays a transition from linear depen- 
dence to a constant value at a critical size. We predict 
that for circularly samples the activation energy will con- 
tinue to increase beyond a critical size, albeit more slowly. 
We have also provided an analytic estimate for the ac- 
tivation energy for nonuniform magnetization reversal in 
nanomagnets with perpendicular anisotropy, which can 
guide experiment and memory device development. 
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Appendix 

The magnetostatic tensor has been calculated by 
Newell [16] . Using Gauss' Law Eq. [3] can be transformed 
into 
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where S is the surface enclosing the sample and n is 
the unit vector normal to the surfaces of the sample. 
The components of the magnetization tensor are ob- 
tained by selecting the appropriate surfaces of orienta- 
tion, e.g. the N xy term is obtained by using the sides 
where ndS = ~kdydz and n'dS' = ydxdz. Notice that 
the trace of the integrand in Eq. [3] can be calculated us- 
ing -V • V'ij^pj) = A(^) = 47n5 3 (r - r'), where A 
denotes the Laplacian. It follows that Tr [N] = 1 for a 
body of arbitrary shape. 
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